Connection between Dispersive Transport and 
Statistics of Extreme Events 



b 



K.W. Kehr a '*, K.P.N. Murthy a - fe , and H. Ambaye a 
a Institut fur Festkorperforschung, Forschungszentrum Jiilich GmbH, 

D-52425 Jiilich, Germany 
Materials Science Division, Indira Gandhi Centre for Atomic Research, 
Kalpakkam 603 102, Tamil Nadu, India 
* E.mail: k.w.kehr@fz-jiilich.de 



Abstract 



A length dependence of the effective mobility in the form of a power law, 
B ~ L 1_ q is observed in dispersive transport in amorphous substances, with 
< a < 1. We deduce this behavior as a simple consequence of the statistical 
theory of extreme events. We derive various quantities related to the largest 
value in samples of n trials, for the exponential and power-law probability 
densities of the individual events. 

I. INTRODUCTION 

Dispersive transport in amorphous materials has been studied for almost three decades 
now, and yet it continues to attract considerable attention. In these studies, charge carriers 
are created at one side of a slab of the material and transported across to the other side. It 
has been realised that the transport phenomena in amorphous media cannot be described 
by standard concepts like uniform drift and diffusional spreading, see for exhaustive 
reviews on this subject and || for a popular account. An important and striking obser- 
vation is the (apparent) dependence of the mobility on the thickness of the material. The 
general features of the mechanism of dispersive transport in amorphous materials are well 
understood; the process is subdiffusive and the delay in the transport occurs mainly due 
to trapping events in localized centers of various energetic depths. A fairly large body of 
theory has been developed to describe several aspects of the transport phenomena; already 
the early theoretical work of Shlesinger and of Scher and Montroll |8j could explain many 
of the experimentally observed phenomena, in particular the length-dependent mobility. In 
this note we shall show that this length dependence of the mobility can be easily understood 
as a simple consequence of the statistics of extreme events. 

The statistical theory of extreme events was formulated for discrete random variables 
more than half a century ago by Fisher and Tippett and Gnedenko ||10|| ; it was extended 
and popularized by Gumbel , . Gumbel described various applications that include 
statistics of extreme floods, droughts and fracture of materials. In the field of condensed 
matter, however, the statistical theory of extreme events has found fewer applications. The 
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principal aim of the theory of extreme events is to obtain the statistics of the largest value 
in a sample of n independent realizations of a random variable. In particular the aim is 
to determine the asymptotic (n — > oo) dependence of the extreme values on the sample 
size n. The random variable in the case of dispersive transport is the residence time in 
the trapping centers. The question then can be posed as follows. Given the distribution 
of trap depths, how does the largest trapping time increase with the number of trapping 
centers; the latter quantity is determined by the thickness of the material, (we shall use 
length synonymously with thickness). Together with the assumption that the extreme events 
determine the behavior of the sum of the residence times, the statistical theory of extreme 
events makes then a prediction concerning the dependence of the largest residence time on 
sample size or the thickness of the material. Since mobility is deduced from the sum of 
the transit times and this quantity is given by the sum of the residence times, the length 
dependence of the mobility would follow in a natural fashion, as we shall show in this paper. 

The paper is organized as follows. In section II the elements of the statistical theory 
of extreme events that are necessary for the later derivations are briefly described. The 
application to dispersive transport in disordered materials is discussed in section III. Detailed 
calculations of various quantities that are relevant for the characterization of the statistics of 
extreme events are made in section IV. The paper closes with concluding remarks in section 
V. 



II. ELEMENTS OF STATISTICS OF EXTREME EVENTS 

This section reviews standard material of the statistics of extreme events [II, 12] which is 



needed later on. Let f(x) be the probability density function (PDF) of a random variable X 
and F(x) = J x dx'f(x') its cumulative distribution function (CDF). F(x) is the probability 
that a particular realization of X has a value < x. Let Q n = {x±,X2, ■ ■ -x n } be a set of n 
independent realizations of X sampled from the density f(x). Let x n i denote the largest 
value of X in the set fl n . The CDF of x n \ is the same as the probability that all the values 
of X in the set Q n are less than x, and hence is given by, 

d> n (x) = F n (x) . (1) 

The probability density of the largest value of X in the set Q n is found by differentiation, 

<P n (x)=nF n -\x)f(x). (2) 

Given a value of u, let g u {v) denote the probability that the first v — 1 values of X are 
less than u and the v- th value is greater than u. An expression for g u {v) can be readily 
written down as 

g u (v)=F»- 1 (u)ll-F(u)}. (3) 
The mean value of v is denoted by n and is given by 

OO 

n={,) = ^g u {,) = Y - Wy (4) 
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The quantity n is the mean number of steps required to exceed a given value of u. The 
above relation can be 'inverted' as follows. For a given n, considered now as a parameter, 
let u n i denote the value of u that obeys Eq. (H). Thus we get an implicit equation for u n i 



F{u nl ) = 1 - - . (5) 



Gumbel fly] calls u n i as the 'expected' largest value of X in a sample of n independent 
realizations. Notice that u n \ is not the mean of the extreme value in a sample of size n. 
Hence u n i should strictly be viewed as a quantity defined by Eq. @. This simple expression 
for the 'expected' largest value helps gain insight into the asymptotic behaviour of the 
extreme values; we shall use this definition of 'expected' largest value to derive the power 
law dependence of the mobility in the next section. 

A simple example is provided by the exponential PDF f(x) = exp(— x) , x > 0. Solution 
of eq. (13) with respect to u n \ yields 

u ni = ln(n) , (6) 

i.e., a logarithmic increase of the expected largest value with sample size. 

There are essentially three classes of behavior of the largest value u n \ for large n, leading 
to different asymptotic forms of $> n (x) for large n [fTT| , p~2|1 . The first class (I) is formed by 
probability densities f{x) which decay at least exponentially for large n. The second class 
(II) results from probability densities whose moments diverge beyond a certain order. The 
third class (III) comprises probability densities where the values of x are bounded. We will 
encounter class I and class II behavior later, depending on the physical quantity that is 
considered. 



III. MOBILITY IN DISPERSIVE TRANSPORT 

In a typical experiment on dispersive transport, a slab of a finite thickness L (usually 
thin films of an amorphous substance) is coated with semitransparent metal electrodes. A 
constant voltage is maintained across the slab. Charge carriers are created at one surface at 
time t = by a laser pulse; the charge carriers are drawn through the slab by the electric 
field. The current I(t) exhibits different behavior at short and at long times, which is most 
clearly seen when plotted on a log-log graph. Shlesinger ]7| and Scher and Montroll f| 
predicted 

U-° t>>l. < 7 > 

where the parameter < a < 1. Deviations from the ideal behavior Eq.(|7|) are still a topic 
of current research (see, e.g. H5,13] and the references therein). The physical explanation of 



the behavior of I(t) is by trapping of the charge carriers in trapping centers with widely 
differing depths. In the short-time regime, most of the charges are within the slab, while at 
longer times the charges are extracted from the slab. 
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A transit time tt r can be deduced from the crossover between the short-time and long- 
time behavior. An effective mobility is then defined by the ratio of effective velocity and 
applied field F, 



t tr F ^ ^ 

In the multiple-trapping model, which is employed here, the transit time tt r is the free transit 
time tf ree plus the sum of all dwell times Tj in the trapping centers. Similar arguments could 
be applied to trap-controlled hopping models. The free transit time is given by the mobility 
B , if there are no trapping events: tf ree = L/{BqF); it is usually short compared to the 
sum of all dwell times. Hence to a good approximation 



' tr 



E 



(9) 



where n is the number of trapping events. For a constant trapping rate the number of 
trapping events is proportional to the thickness L of the slab. 

For broad distributions of dwell times Tj, the sum in Eq.(|9|) should be dominated by the 
largest dwell time. With this assumption 

Ur ~ T largest ■ (10) 

The Arrhenius law is assumed for thermally activated processes, 

n = r exp(^-^) (11) 

where E{ is the energy required for release from the trapping center i. The largest dwell time 
is then determined by the largest trapping energy. The probability density of the depths of 
the trapping centers, i.e., of the energies necessary for release is assumed to be exponential, 

fXE)=E- 1 e* V (-^ E>0. (12) 

Note that f(E) is normalized to unity in the interval (0, oo). It is easy to convert the 
probability density of the trapping energies into the probability density of the dwell times 
p(r), using the Arrhenius law Eq . fllTD . We have 

p (r) = T <T<OO (13) 

with the parameter a = kBT/E c . The probability density (|13"D is normalized, but already 
its first moment does not exist for < a < 1, which is the parameter range of interest for 
dispersive transport. The cumulative distribution function of r is given by 

P(r) = 1 - (-)- • (14) 

To 

Equation (|5]) for the expected largest value in a sample of n trials yields 
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1_ 

Tnl = T n<* 



(15) 



As already stated, the average number of trapping events n is proportional to the thickness 
of the experimental sample L. Hence we expect, using the assumption that the transit time 
is dominated by the largest dwell time 

t tr ~ L« 

B ~ L 1 -^ . (16) 

This is precisely the behavior of the mobility that has been observed in experiments on 
dispersive transport, see for instance ||. Here we have derived this behavior from the 
statistical theory of extreme events. 

IV. EXPONENTIAL AND POWER-LAW PROBABILITY DENSITIES 

A. Motivation 

Various questions arise with regard to the validity and the significance of the above result 
for the length dependence of the mobility. For instance, the meaning of the "expected largest 
value in a sample of n trials" is partially intuitive; i.e. the precise meaning of the quantity 
that follows from Eq.(^) is different from what is suggested by this notion. Precisely defined 
quantities are the moments of the probability density function of the largest value in a sample 
of n trials (if they exist) and the most probable value. Of course, complete information is 
contained in the PDF itself. Fortunately, all quantities of interest can be derived exactly, 
if the underlying PDF for one event is exponential, or of power-law form. This section will 
describe the results of these calculations, including a numerical determination of the PDF 
of the dwell time. 

B. Exponential probability density for single events 

The basic quantity which determines the dwell times of particles in the multiple-trapping 
model for dispersive transport is the trapping energy E. It is a random quantity and the 
simplest, experimentally relevant, assumption is the exponential PDF, see Eq. (|12|) . The 
dimensionless form of this PDF is f(x) = exp(— x) with the variable x = E/Ec and the 
restriction < x < oo. The cumulative distribution function is then F(x) = 1 — exp(— x). 
The expected largest value in samples of n trials has already been given in eq.(|6|). The PDF 
for the largest value x n i in a sample of n trials follows from Eq.(^) as 

^\x)=n(l-e-*) n - 1 e- x . (17) 

The index (1) shall indicate that this PDF is of type I in the classification of Gumbel |ll| , [l2| . 

The moments of (p^(x) can be calculated exactly. The first moment, or mean value of 
x n i is defined as 

POO 

(x n i) =n dxx(l- e~ x ) n - x e~ x . (18) 

J 
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The evaluation of the integral is made in the Appendix. The result is 

(x n ) = V(n + 1) - • 
The digamma function ip{z) is defined as [14| 

where 7^ is the Euler constant. For integer n 

ip(n+ 1) 



ip(z) = — In (T(z)) and 



-IE 



j=i J 



Thus 



n 1 

3=1 J 



(19) 



(20) 



(21) 



(22) 



Euler's constant is asymptotically given by 



IE 



lim 
n— >oo 



n 1 
3=1 J 



ln(n) 



Hence the first moment is given for large n by 

(x n l) =1e + ln(n) 



(23) 



(24) 



Note that the first moment differs from the expected largest value of x n i by Euler's constant, 
cf. Eq.(^|). Figure 1 contains the results on the first moment Eq.(|22|) and the expected 
largest value Eq.(||). Results on a numerical determination of the mean value (x n i) have 
been included, to demonstrate that this quantity can also be accurately determined by 
numerical simulation. 

The second moment (x^i) can be calculated by the same technique as employed for the 
first moment, cf. the Appendix. The result for the second moment is 



(x 2 nl ) = mn + 1) - ^(l)f + - 4> {1 \n + 1) 

where ip^ l \z) is the threegamma function defined by [14 



d 



^{z) = ^- 2 {\nT{z)) . 

The variance of the extreme value is then given by 

(x 2 nl )-(x nl ) 2 = ^(l)-^(n + l). 
Therefore, the variance of the extreme value is 



a 



n 1 00 1 

E F = «2)- E u- 



k=l 



k=n+l 



k 2 



(25) 



(26) 



(27) 



(28) 
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Note that £(2) = tt 2 /6. Hence the variance of the extreme value approaches asymptotically 
this value, i.e., it is asymptocially independent of n. The significance of this fact has been 
stressed by Gumbel [O]. It is easy to derive the following bounds for the variance: 

7T 2 1 2 7T 2 1 

<a 2 < . (29) 

6 n 6 n + 1 

The variance of the dimensionless extreme energy as a function of n has been included in 
Fig.l. 

We return to the PDF for the largest value x n \ in a sample of n trials, which is generally 
given by Eq.(|2]) and for an exponential PDF of single events by Eq. flTTD . We define the 
most probable value of the random variable x n i , i.e., the most probable extreme value, as 
the value of x where the PDF tp n (x) has its maximum. Differentiation of Eq.(0) yields the 
extreme value condition 

<p' n (x) = n {F n -\x)f{x) + (n - l)F n -\x)f{x)f{x)} = . (30) 

This equation can be cast into a simple form, 

For an exponential PDF f(x) = exp(— x) the solution of this equation yields the most 
probable extreme value x m _ p _ as a function of the sample size, and is given by 

x m . p . = ln(n) . (32) 

Thus the most probable extreme value agrees with the expected extreme value in the case 
of the exponential PDF for the single events, cf. Eq.@. 



C. Power-law probability density for single events 

The quantity that determines the effective mobility is the sum of all dwell times, which 
is assumed to be dominated by the largest dwell time. In this subsection we direct the 
attention to the statistics of the largest dwell time. The dwell time in a trapping center 
follows from the trapping energy via the Arrhenius relation Eq. (|i~l"D . The PDF and the CDF 
for the single events were already given in Eq.(|l3"D and Eq. (|T4]) , respectively. We will use 
dimensionless variables y = t/t henceforth. The PDF for the extreme value y in a sample 
of n events is then given by 

^\y) = nail - y~")^y-^ 1 < y < oo (33) 

where the parameter a = IzbT / E c and the relevant range is < a < 1. The index (2) shall 



also indicate that the distribution of y is of type II in the sense of Gumbel |11]JT! 

The moments of ^\y) can be calculated. Let (y^) denote the kth moment of the y n \. 
It is easily shown, 

(y k nl ) = n(3(l--,n) (34) 
a 
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where f3(a, b) is the usual beta function. 

^From the properties of the beta function it is clear that the kth moment of the extreme 
value exists only if k is less than a. Since in our problem, a is between zero and one, no 
integer moment k > 1 of the extreme value exists. Nevertheless, the expected extreme value 
(as defined by Gumbel) exists, from which we have derived the length dependent mobility. 

It is instructive to determine the most probable extreme value from the PDF Eq. (|33| ) for 
the largest value y n i in samples of n trials. The condition for an extremum has been given 
in Eq.(^); evaluation with the PDF for single events Eq. (|T3|) and the CDF Eq.(|14|) gives 

Vm.p. = ~z ' • (35) 



1 + a 

For large n, we can replace 1 + an by an, and we get, 

^. = (t^) 1/ °" 1, °- < 36 > 

Another quantity, which is often used to characterize broad probability densities, is the 
typical value. It is defined by 

hi(x typ ) =< ln(x) > . (37) 

where the brackets indicate a mean value taken with a general PDF. The typical value of 
the PDF Eq.(|33|) for the extreme value of r is given by 



dyhi(y)naO--y~ a ) n ~ 1 v~ a ~ 1 - (38) 

By the substitution z = y~ a this integral is transformed into 



71 



i 

n-l 



\n(y typ ) = -- J q dz\n(z)(l - z)^ 1 (39) 

The above integral can be evaluated exactly employing essentially the technique described 
in the appendix (for the exponential PDF) and the result is 

ln(^ P ) = --[^(n + l)-^(l)]. (40) 

a 

Asymptotically, the typical extreme value behaves as 

y tyP — >n*exp(— ). (41) 
a 

That is, y typ is related to < x n i > by y typ = exp(-< x n i >). We point out that the connection 
of the typical extreme value of y to the mean extreme value of x is generally valid, if these 
two variables x and y are related by the exponential transformation y = exp(x). Notice 
that the energy E and the dwell time r are related to each other by the Arrhenius law (11), 
which is an exponential transformation. 

In Fig. 2 we have plotted the expected largest value r ni /r of the dimensionless dwell time, 
the most probable value according to eq.fllOD, and the typical value following from Eq.(E0), 
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for a = 0.5. In the case of the dwell time, the expected, the most probable, and the typical 
largest value are all proportional to n~ , with different n-independent factors for large n. It 
seems to be a general fact that the quantities which characterize the largest value in samples 
of n trials, show similar behavior with respect to n, if they exist. 

Although the PDF for the extreme value of y in samples of n trials is exactly given 
by Eq.(^), the evaluation of this function is not practical for large n. Therefore, we have 
determined this PDF for fixed values of n by numerical simulations. We have generated 
random energies according to the distribution Eq.(^2|) with E c — 1 and calculated dwell 
times using the Arrhenius law (|TT|) with r = 1 and a = 0.5. The normalized dwell times 
were sorted into bins containing L values. The inverse of the lengths of the bins gives the 
PDF, when properly normalized. The result for n = 1024 is given in Fig.3. The maximum 
of the distribution agrees with the prediction Eq. fl35|) . The typical value is found right of 
the maximum, close but not identical to the median value. Similar observations have been 
made for broad distributions previously, for instance in Ref. P"3 . 



Note that the PDF of the extreme value can have a different form, if it is determined by 
different techniques. If the dwell times are binned into intervals of logarithmically increasing 
intervals, the maximum of the resulting PDF is found at a different location. This is due 
to the fact that by this procedure another PDF is estimated, which is related to the one we 
used by the usual transformation with a Jacobian. The maximum of the latter distribution 
is at ria which is identical to the expected largest value as introduced by Gumbel. The 
message is that the precise characterization of an extreme value depends on the PDF used 
for this purpose. 

A final problem is the justification of the replacement of the dwell times, cf. Eq.(|9|) by the 
largest dwell time, as was done in Eq. flToD . In the subsequent derivation, the largest dwell 
time was identified with the expected largest time, however, any quantitiy that characterizes 
an extreme value can as well be used in the argument. The replacement can be justified by 
extending the derivations given above to the most probable K th extreme value. Let v 7 n,^(^) 
denote the probability density function of the K th extreme value. A formal expression for 
(fin^ix) can be easily derived and is given below: 



^) = { n-K)l(K-l)l f{x)Fn ~ K{x)[1 ~ F{X)]K " (42) 

Note that if we set K = 1, we recover the probability density of the first extreme, which we 
have considered so far. The condition for an extremal value is 

<p' n>K (x) = . (43) 

Solution of this equation with respect to x gives the most probable K th extreme value. 

The extremum condition Eq.(|4"3|) can be solved for a power- law PDF as given in Eq. (|T3|) . 
The result for the dimensionless variable y = t/tq is 



1 + an 



« = lraJ • (44) 

Let tk denote the ratio of the most probable K extreme to that of the most probable 
(first) extreme. It is given by 
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(45) 



For a = 0.5 and K = 2, we have r = 9/16. The ratio becomes rapidly smaller for larger K; 
it is already small for K = 2 when a is small. The consequence is that the inclusion of the 
second, third, etc. extreme value in the estimate of the dwell time would not change the 
asymptotic dependence on the sample size n; only numerical factors would be modified. 



V. CONCLUDING REMARKS 

The apparent anomalies in the transport properties of charges in amorphous substances, 
which are summarized as 'dispersive transport', are due to broad distributions of trapping 
times of the charge carriers. It is satisfying that features like the length dependence of the 
effective mobility can already be derived from the statistics of extreme events, the extreme 
events being the occurrence of particularly long trapping times. 

Although the length dependence of the effective mobility could already be deduced from 
the notion of the "expected largest value" of the trapping time, we investigated in detail 
various quantities in the context of the statistics of extreme events. Explicit expressions 
could be derived for the mean, the most probable, and the typical extreme values in the case 
of exponential or power-law probability densities for single events. We could also justify 
the use of the largest dwell time to estimate the summary transit time by considering the 
second, third, etc. extreme values. Their effect is to modify numerical factors, but they do 
not alter the asymptotic dependence on the sample size. 

All derivations were made here for the ideal case of an exponential density of energy 
levels. In practice deviations from the exponential density of states are of interest. It turned 
out that many features of dispersive transport are also present for a Gaussian density of 
states ||; if the temperature is sufficiently low. Hence it is of interest to extend the present 
derivations also to the case of other densities of states, for instance to the Gaussian one. In 
the case of more complicated probability densites of the single events, analytical derivations 
are either difficult or impossible. Hence one has to resort to numerical simulation to study 
the statistics of extreme events in those cases. 

Finally we point out that the observation of different behaviors of distributions, depend- 
ing on which of two exponentially related variables are used, is rather general. For instance, 



de Arcangelis, Redner, and Coniglio [16] found broad distributions of voltage drops in ran- 
dom resistor networks, whose moments are characterized by an infinite set of exponents. To 
the contrary, the distribution of the logarithm of the voltage drops behaves normally, in that 
the moments exhibit constant-gap scaling. Similar behavior of the probability distributions 
of the cluster numbers in the percolation problem was discussed by Stauffer and Coniglio 
11711 - 



Discussions with G. Schiitz are gratefully acknowledged. KPNM is thankful to M.C. 
Valsakumar for numerous discussions on the statistics of extreme events. 
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VI. APPENDIX: EVALUATION OF FIRST MOMENT OF <p$(X) 

In the expression for the first moment Eg . (|i"8|) the substitution y = exp(— x) is made, 

(u») = - tdynlniy^l-yT' 1 . (46) 



ivm y e - 1 



The following trick |18| is applied: 

Hv) = e^O^ (47) 
We have 

(«n> = e^o - n £ dy - y)^ 1 . (48) 

The integral can now be performed with the result 

(u n ) = - n ^ l + ^-^ n) (49) 



where (3(a,b) = Jq 1 dx x a 1 (1 — x) b l . Hence we have 



(u n ) = -n^-(3(z,n)\ z=1 . (50) 



The beta function is related to the gamma function (see |T 



0M) = ™ (M) 
Differentiating and observing that T(n + 1) = nT(n) and T(l) = 1 we obtain 

{Un) ~T(n + l) T(l) ' {b2) 
Using the definition of the digamma function |I4J this can be written as 

(u n ) = £(n + !)-£(!)■ (53) 



11 



REFERENCES 



[1] G. Pfister and H. Scher, Adv.Phys. 27, 747 (1978) 

[2] E.W. Montroll and B.J. West, in Fluctuation Phenomena, edited by E.W. Montroll and 

J.L. Lebowitz (North Holland, Amsterdam, 1979), Chap. 2. 
[3] F.J. Schmidlin, Philos.Mag.B 41, 535 (1980) 

[4] T. Tiedje, in The Physics of Hydrogenated Amorphous Silicon, edited by J.D. Joanopou- 
los and G. Lucovsky, Topics in Applied Physics Vol. 56 (Springer, Berlin, 1984), Chap. 
6. 

[5] H. Baessler, phys.stat.sol.(b) 175, 15 (1993) 

[6] H. Scher, M.F. Shlesinger, and J.T. Bendler, Physics Today, January 1991, p. 26 
[7] M.F. Shlesinger, J.Stat. Phys. 10, 421 (1974) 
[8] H. Scher and E.W. Montroll, Phys.Rev.B 12 2455 (1975) 
[9] R.A. Fisher and L.H.C. Tippett, Proc.Cambridge Phil.Soc. 24, 180 (1928) 
[10] B.V. Gnedenko, Ann.Math. 44, 423 (1943) 

[11] E.J. Gumbel, Statistical Theory of Extreme Values and Some Practical Applications, 

National Bureau of Standards Applied Mathematics Series 33 (1954) 
[12] E.J. Gumbel, Statistics of Extremes, (Columbia Univ. Press, New York N.Y., 1958) 
[13] A. Jakobs and K.W. Kehr, Phys.Rev.B 48 8780 (1993) 

[14] Handbook of Mathematical Functions, edited by M. Abramowitz and LA Stegun, (US 

Government Printing Office, Washington D.C., 1965) 
[15] K.W. Kehr and K.P.N. Murthy, Phys.Rev.A 41, 5728 (1990) 
[16] L. de Arcangelis, S. Redner, and A. Coniglio, Phys.Rev.B 31, 4725 (1985) 
[17] D. Stauffer and A. Coniglio, Physica 143A, 326 (1987) 
[18] M.C. Valsakumar, unpublished notes (1977). 



12 



Figure Captions 



Fig. 1: Mean and variance of the dimensionless largest energy in samples of size n, taken 
from an exponential PDF, as functions of n. Full line: expected largest value x n i, dotted 
line: mean of largest value, dash-dotted line: variance of largest value. Symbols: O : mean 
largest value by simulations, + : variance by simulations. 

Fig. 2: Normalized extreme dwell times as functions of the sample size n. Full line: most 
probable largest dwell time, short dashes: expected largest dwell time, long dashes: typical 
largest dwell time. 

Fig. 3: Probability distribution function of the largest value of the normalized dwell time 
in samples of size n = 1024. The PDF was determined from 10 5 realizations. Indicated are: 
most probale larges value (full line), expected largest value (short dashes), median (points), 
and typical largest value (long dashes). 
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